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Abstract 



The optimal (Monge-Kantorovich) transportation problem is discussed from several 
points of view. The Lagrangian formulation extends the action of the Lagrangian L(v, x, t) 
from the set of orbits in R" to a set of measure- valued orbits. The Eulerian, dual for- 
mulation leads an optimization problem on the set of sub-solutions of the corresponding 
Hamilton- Jacobi equation. Finally, the Monge problem and its Kantorovich relaxation 
are obtained by reducing the optimization problem to the set of measure preserving map- 
pings and two point distribution measures subjected to an appropriately defined cost 
function. 

In this paper we concentrate on mechanical Lagrangians L — \v\ 2 /2 + P(x, t) leading, 
in general, to a non-homogeneous cost function. The main results yield existence of a 
unique flow of homomorphisms which transport the optimal measure valued orbit of the 
extended Lagrangian, as well as the existence of an optimal solution to the dual Euler 
problem and its relation to the Monge- and Kantorovich formulations. 

1 Introduction 

1.1 Motivation 

Consider the Schrodinger equation on a domain Q C M n 



where P = P(x,t) is a given potential (pressure field). If d£l ^ we specify a boundary 
condition under which equation Q1.1J1 is well posed and generate a unitary semigroup. We are 
interested in the following question: Given a pair of probability measures po,p\ on Q (say, 
absolutly continuous subjected to the densities po, pi), can one solve the two point boundary 
problem, given only the data po(x) = \ip\ 2 (x,0) and pi(x) = \i/j\ 2 (x,T)? 

This problem is a prototype for a variety of wave equations admitting a Lagrangian 
formulation. However, in spite of the underlying Lagrangian stracture associated with 
we do not yet know wether such a two point boundary value problem is well posed for this 
equation. 

We turn, therefore, to the semi- classical approximation h — > 0, where ip is written as 
^fpe % ^l h . Then the Schrodinger equation is reduced to the Hamilton-Jacobi (HJ) equation 
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coupled with the continuity equation 

Pt + V x ■ [pVx<t>] = • (1.3) 

The two point boundary problem for the Schrodinger equation in the semi classical approxi- 
mation is, therefore, reduced to solving (|1.2|) and Q1.3JI under the sole data of the initial and 
end probability densities po = p{,0) and p\ = p(,T), without any data on the velocity field 
V x whatsoever. 

We attempt to associate a Lagrangian on the state space of orbits of probability densities 
p(,t) defined for each t £ [0, T]. Since we wish to generalize the standard Lagrangian for- 
malism defined on classical orbits x : [0, T] — > £1, we shall extend this to a space of "relaxed 
orbits" of the from [x = fj,^(dx)dt, where p^t) is a Borel probability measure on £1 for any 
t 6 [0,T]. The special class of deterministic orbits is embedded in this set as p^ = 5^^), 
where x is a classical orbit. 

The next step is to define a metric on this set of relaxed orbits which generalizes the 
action \x(t)\ 2 dt for deterministic orbits. It turns out that the above problem is closely 
related to the classical Monge-Kantorovich problem of optimal transportation subjected to 
quadratic cost. An historical background of this is given below. 



1.2 Historical Background 

The classical problem of optimal mass transportation was suggested by Monge in the 18'th 
century [M]: given a cost function c{x, y) (originally, c = \x— y\) and a pair of Borel probability 
measures po, p\ on (say) a common probability space $7, minimize 

J c(x,T(x))vo(dx) (M) 

along all Borel mappings T : 0, — > which transport po into p± (T#/xo = /xi), namely 

po (T _1 A) = pi(A) V Borel sets A C £1 . (1.4) 

The Monge problem was revived in the last century. In particular, Kantorovich [K] introduced 
in 1942 a relaxation, reducing the Monge problem to a linear programming in a cone of two- 
point distributions over Q whose marginals are po,Pi respectively: 



mm 

A 



c{x,y)\(dx,dy) ; n$ A = p , n^'X = pi (K) 



Here Tr^,i = 0,1 are the natural projections of Q x Q on its factors. A particular attention 
is given to the Wasserstein metrics 



W p {no,p,{) 
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A 



\x - y\ p X(dx,dy) ; n£ X = p , tt# A = pi 



i/p 



(1.5) 



where p > 1. 

The most striking advantage of the relaxed problem is that a minimizer always exists 
by the compactness of the set of probability measures (assuming c is continuous and f2 is 
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compact). It can be shown [Am] that, if c is continuous and contains no atoms, then 
the minimum of the Kantorovich problem coincides with the infimum of the Monge problem. 
The existence of an optimal Monge mapping is reduced to existence of such a minimizer of 
the Kantorovich problem which is supported on a graph of a Borel map. 

Further progress was achieved in the last few decades. In the late 80's and early 90's 
Brenier [B] studied the Monge problem with a quadratic cost c = \x — y\ 2 on R n and showed 
the existence of a unique, optimal Monge map provided Ho,Hi posses finite second moments 
and fiQ is absolutely continuous with respect to Lebesgue. Moreover, he showed that this 
map is the gradient of a convex function <& which solves the dual problem 

infj f $/x (dx)+ f Vdfii(dx)\ ; <l),$eC(l") ; + > x ■ y V x, y G R n , 

(1.6) 

It turned out that any map T which transports hq to fi\ must be of the form V$ o S 
where S perseveres hq. This results is interpreted as a Polar factorization for mappings, 
generalizing the matrix polar factorization. Another interpretation of V*!* is as a monotone 
vectorized rearrangement in the class of maps transporting to fi\. A generalization of this 
result also holds for more general, strictly convex (and concave) homogeneous cost functions 
c(x,y) = h(x — y), using special definitions of convexity ([GM], [C] and references therein). 

1.3 Objectives and main results 

In general, if fio contains an atom, then there is, in general, no deterministic mapping T of 
any type which maps no into /xi, so there is no sense to compare the deterministic Monge 
problem (M) with the probabilistic Kantorovich problem (K). However, we may still consider 
the following alternative formulation in terms of an optimal flow with respect to some family 
of cost functions ct lt t 2 = J(x,y,ti,t2)' 

(F) : Find a relaxed orbit fi = H(t)dt and a flow of diffeomorphysms : — > f2 for 
h,t 2 e (0,T) such that 

(i) is the optimal Monge mapping with respect to Ct ly t 2 transporting to fj,^ for 
anyh,t 2 G (0,T). 

(ii) lim t ^o = fj,Q and lim t _>T = Hiin the weak sense of measures. 

(Hi) The limits lim^y T* =: exists uniformly and T^is a continuous mappings for any 
rG(0,T). 

It is feasible that, once a solution to the flow problem F is provided, a co,r optimal 
solution to the Monge problem M with respect to no, n\ exists by T = lim T _>o provided 
the later limit exists as a Borel map. 

Our starting point is the definition of a norm of a measure- valued orbit as the min- 
imal L^-norm of the velocity fields v which satisfy the weak form of the continuity equation 

j« = «(x,t); £ Jjci>t + v-V x (t>]L H t)(dx)dt = ; V</> £ Cq(£1 x [O.T])| . (1.7) 
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and 



inf / \v\ P H( t \(dx)dt 
v Jnx[o,T] 



i/p 



where the infimum is taken over all //—measurable vectorfield v satisfying (|1.7|) . Denote 
the set for which \\fi\\ p < oo as H p . This is a normed cone. In section |2] we shall indicate 
some of its properties and prove a compactness embedding of H p (for p > 1) in a set of 
orbits which satisfies Holder continuity in an appropriate topology. In particular, the end 
conditions fio := M(o)> Mi := M(t) are uniquely defined for fi £ H p where p > 1. 

In the rest of the paper we concentrate on the case p = 2. The connection between the 
cost function ct 1} t 2 = J posted in formulation (F) above and the pressure P is as follows: 
The function J = Jp is the action associated with the Lagrangian 



J P (x,y,ti,t 2 ) = inf 



1-2 



\m\ 2 



+ P(x(t),t) 



dt ; x : [ti, £2] — ► ^7 x(ti) = x , £(£2) = V 



The main result of this paper, formulated in section |3J reveals a connection between the 
following approaches: 

L: The Lagrangian approach: Minimize a Lagrangian Lp on the space of orbits /i G H2: 
£ := inf L P (/x) ; £p(//) := ^|Ml2+ [ / P(i(t)(dx)dt , fi £ H2 , M(o) = Mo, M(T) = Ml 

^ JO J!l 

E: The Eulerian approach: Maximize on the set of velocity potentials <j) 



£ := sup 



4>(x,T)m(dx) - / (j)(x,0)no(dx) 
n in 



where the supremum is taken in the set of all functions <p = <fr(x, £) which are sub- 
solutions of the Hamilton-Jacobi (HJ) equation (|1.2j) in a sense to be defined. 

M: The Monge approach: Minimize on the set of mappings verifying (|l,4j) 
M :=inf I J J P (x,T(x),0,T)fi (dx) ; T #/ u = Ml} 

K: The Kantorovich approach: Minimize on the set of 2-point probability measures with 
prescribed marginal 

/C:=min|y J J P (x,y,0,T)X(dx,dy) ; ir^X = /x , n^X = Mi} 



Our first result reveals the relation between the above formulation: If P £ C l (Q. x [0,T]) 
then 

£ = £ = IC 
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holds for arbitrary (probability, Borel) end measures hq,/j,i. As discussed above, the Monge 
problem may not have a solution at all (e.g., if /xo contains an atomic measure and the set of 
transporting mappings T#/i = [i\ is empty). 

The second part of our main result shows the relation between the flow problem (F) 
and the Lagrangian formulation L. This is the relation between the optimal velocity field v 
realizing ()1.8|) and the induced flow 

jT t tl (x)=v(T t ti ,t) . (1.10) 

To elaborate, we shall prove 

1) There exists a minimizer [i £ H2 of L which satisfies the end conditions. This minimizer 

may be non-unique. 

2) There exists a maximizer ip of E which is a Lipschitz function on ft x [0, T] and satisfies 

the equation 

i> t + \\v x M 2 = P (i-ii) 

almost everywhere. Again, such a maximizer may be non-unique. 

3) The vector field v = V^V is defined everywhere on some relatively closed set Kq C 
ft x (0,T) which contains the support of any minimal path /i of L given by (1). 

Under some additional assumption on P (see Main Theorem in section |3J) we also get 

4) The vector field v = is locally Lipschitz continuous on Kq. 

5) The restriction of v to the support of any minimal orbit of L is uniquely determined. 

6) The flow T induced by v Q1.10|l leaves Kq invariant. 

7) The flow transports to //( t2 ) for any minimizer \x of L and any t\,t2 € (0, T). 

Moreover, it is an optimal Monge transport with respect to the action Jp(-, ■,t\,t2). 

8) The maps lim r _ T T[ := Tf : £1 -> Q and lim^oT* := T Q : ft -> Q exist and are 
continuous for any t £ (0,T). Moreover, [Tj 1 ]^ (res. [Tq]_^) is an optimal Monge map 
with respect to the action Jp(,,t,T) (res. Jp(-,-,0, t)) transporting fi^ to fi\ (res. hq 

to /i (t) ). 

9) If limt_> T Tq := T exists Borel map, then T transports fiQ to m and is an optimal 

solution of the Monge problem M. In this case 

M = C = £ = K 



A particular case is the pressureless flow P = 0. Here the optimal potential satisfies 

^ + ^|V^| 2 = (1.12) 
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and the associated action is 

I 12 

\x — y\ 

Jo{x,y,ti,t 2 ) = Try- -r , 

reducing the Monge-Kantorovich problem to the Wasserstein metric W2 for quadratic costs 
(|1.5|) . The associated flow, claimed in (6), is given in this case by 

T t t \{x) = x + (t 2 -t 1 )V x ^(x,t l ) 

where V x ip is defined and Lipschitz everywhere. In particular it follows that, for a quadratic 
cost, an optimal Monge map T#/xo = exists and is unique provided V x ip(x,0) is 
measurable. 2 In this case, Brenier representation T = = V x <& of the optimal map [see 
(|1.6I) and the proceeding discussion] is recovered via 

<Z>(x) =x 2 /2 + T^(x,0) . 

The connection between the Monge-Kantorovich problem in the quadratic case and the flow 
problem L (P = 0), as well as the dual relation E together with the Hamilton- Jacobi equation 
Q1.12J1 was indicated by several authors (see [BB], [BBG]) 3 as well as in the excellent mono- 
graph of Villani [V]. However, to the best of my knowledge, the existence and uniqueness 
result for the flow without any regularity assumptions on the end measures fio,H\ is new 
even in the case P = 0. In fact, the existence and uniqueness of the flow holds even if there 
is no optimal Monge map. 

In section |3] we shall start to develop the tools needed for the proof of our main results. 
Section 14.11 deals with a dual formulation for the norm \\[i\\2 for an orbit of measure [x = 
fimdt £ H2. It follows that 



J J\V x <f>\ 2 fi(dxdt) 



' sup 



where the supremum is taken on the set of test functions <fi(x,t) = <j) € Cg(S7 x [0,T]). An 
equivalent definition turns out to be 



1 



|=sup^- / P(x,t)n(dxdt) - / (j)(x,0)fi (dx) + / 4>(x,T)n\(dx) \ (1.13) 
4>,p I Jfti Jn in 



where the infimum above is on the pairs of "velocity potentials" <j) e C 1 ^ x [0,T]) and 
"pressures" P = P(x.t) which are related via the Bernulli-type (or Hamilton- Jacobi) equation 
()1.2|) . In case of a prescribed pressure P (as in this paper), this identity reveals the relation 
between the Lagrangian formulation L and the Eulerian one E. In section l4~2l we imply a dual 
formulation to a strict convex perturbation of the Lagrangian L, leading to an approximation 
of the Euler formulation E, to be used in the proof of the main result. 

For the proof of the main result we shall also need a series of auxiliary Lemmas and 
definitions related to the Hamilton-Jacobi equation. In subsection 15 1 1 1 we list these definitions 

2 Since ip(,0) is a Lipschitz function, V x ip(x, 0) is a measurable function defined a.e, so we recover the 
existence of an optimal map if /10 is a continuous w.r to Lebesgue measure. 

3 I wish to thank Prof. D. Kinderlehrer for turning my attention to these publications. 
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and Lemmas, concerning forward (maximal), backward (minimal) and reversible solutions of 
the Hamilton-Jacobi equation, which are essential to the proof of the main results. The proofs 
of the Lemmas are given in subsection 15.21 In 15.31 we utilize these results for the proofs of 
our main Theorem. 

In the rest of the paper we shall restrict ourselves to the flat torus f2 = IR n /Z n . The 
reason is that we wish to avoid compactness problems originated from measures on M n , on 
the one hand, and the boundary conditions for the Hamilton-Jacobi equation required in 
case of a bounded domain f2 C W n . The flat torus is the simplest example in the sense that 
it is compact manifold with no boundary, on the one hand, and it inherits the Euclidean 
geometry from R n on the other. Any function (or probability measure) on is understood 
as a periodic function (or periodic, normalized per-period measure) on R n , unless otherwise 
is explicitly specified. In particular, a mapping T : £1 — * £1 is understood as a mapping on the 
covering M. n which satisfies T(x + z) = T(x) + z for any x £ R n and any zeZ", 
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List of symbols and definitions 

• n := R n /Z n . 

• J =[0,21 ; / = (0,T) 

• Jlj = Q x 2 , £li = Q x 2o. 

• L2P; is the set of all locally Lipschitz functions in f2)2 . 

• M is the set of all probability Borel measures supported in fi. 

• Mi is the set of all Borel probability measures supported on S7/ which are decomposable 
as \i G -4=>- \i = /U( t )dt where G a.e. tel. 

• if /u is Lebesgue continuous measure, then G L 1 (r2/) is the density of fi. 

• 7r^ (res. Tr^ 1 ) ) is the natural projection of Q x $7 on its first (res. second) factor fi. 

• For any pair /j.q, fi\ G Al, the Wasserstein-p metric is defined by 

Wp(fj, ,m) := inf / / |x - y\ p \(dxdy) 
where the infimum is on all probability measures on f! x H such that tt^A = MOj 

(1) X 
TT y # '\ = jJL l . 

• E M (V>) := Jq J n tp{x,t)n {t) {dx)dt. Likewise, E M(t) (V') = f n ii(x,t)n {t) (dx). 

• A lifting v of fi £ A4j is a Borel measure on 0/ x M n such that 

: = / / / V>(^ t)u{dxdtdv) = Epty) ; E„(Vt + v • V^V) = 
for all $ G Cq 1 ^/). 



2 A metric space for measure's orbits 

We start with the following 

Definition 2.1. 2ei fi G A1/. T/ien /i G H p (I,A4) if there exists a lifting v of \i such that 
¥, u (\v\' p ) < oo. We shall also define the H p norm of fx G H p by: 

\\n\\ p =mf[EMn} 1/p 

where the infimum is taken over all liftings of fj,. 
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Lemma 2.1. H p is complete and locally compact under the weak C* topology if p > 1. 
That is, for any bounded sequence [i n in H p we can extract a subsequence which converges in 
C*(Qj) to some \i G H p . In addition: 



lim \\[i n \\ v > 

n—>oo 



Proof. By definition there exists a set of liftings u n corresponding to fj, n . Moreover, this 
sequence can be chosen so that E^(|u| p ) < C, so u n and vi> n are tight on Qi x R n (since 
p > 1 and Qj is compact). Hence the weak limit v of v n is a lifting of the weak limit //of fj, n , 
and E z/ (|v| p ) < C, hence [i G H p . The same argument also yields the lower-semi-continuity 
of Hp. □ 

Lemma 2.2. If \i = fimdt G H p; p > 1 i/ien £/te mop t — * fiu\ is a Holder (p — l)/p 
continuous function from I into Ai with respect to the weak (C* ) topology equipped with the 
Wasserstein-1 norm W\: 



Wi(// 0) A*i)= SU P / <t>{m{dx) - no(dx)) . (2.1) 
|v<£|<i Jn 



Proof. We know that an optimal lifting v exists for \i G H p . The measure v can be decom- 
posed, by the Theorem of measure's decomposition [AFP], into v = fiu\(dx)v Xt t(dv)dt, for \i 
a. a. (x,t). We may define now the velocity field 

v(x,t) = E v (v) 



for fj, a.a. (x,t). It follows that v Gh^ and, moreover, 



By assumption: 



/ Jn 



Ot 



V-(t)(dx)dt 



\v\ p fj,(dxdt) 



i/p 



i Jn 



v ■ V ' x (j)iit t )(dx)dt 



(2.2) 



where 4> = <j)(x,t) is in C£(Q/). Let </>(x,t) = h(t)$(x) with $ G C 1 ^) and 7i G C£(I). Then 
:= J^^(x)n^(dx) satisfies 



•70 



/i(i) / V$ ■ Vfj,($(dx)dt 



n 



By Holder inequality 



/" h(t) ( V$- Vfir t )(dx)dt < |V$| 
o in 



/i 9 (i)(ft 



1/9 



I /'''I 



with g = p/{p— 1). It follows that /$ G W 1,P (I) and, moreover, ||/$||i iP < C| |<3?| |i j00 . This 
implies the result by Sobolev imbedding together with the dual formulation of the W\ norm 

G3J. □ 
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Given /xq and p\ G M, define the set 

A p (/io,Mi) : = = M(t)<& G Hp ;/i (0 )=/i , M(T)=Mi;} • 

Corollary 2.1. T/ie sei A p (/io, where p > 1 is closed and locally compact in C(I; C*(f2)). 

Similar versions of the Lemma and Proposition below can be found in [Am]. We also note 
that Proposition 12 . 1 1 in the case p = 2 is a special case of our main Theorem (see section |3J). 

Lemma 2.3. ( Regularization Lemma): If fi G H p then there exists a sequence \x e G H p 
of smooth density so that p = lim^o M s holds in C*(fij) and, moreover, 



lim I \a £ \ L = I lul 

e-+0 M p 



In addition, for any to,t\ G /, 



limWp^,^) = W p {n {to) ,n {tl) ) . 

We next consider the relation between H p and the optimal solution of the Kantorovich 
problem. 

Proposition 2.1. Assume p > 1. Let no,fJ-i G M.. Then A p (/io,/ii) ^ 0. and 

inf ||/i||p = Wp(// , • 
MGA p (/io,A«i) 

The proof is similar to the proof of Theorem 4.2 of Ambrosio [Am] for the metric case 

09=1). 

We note that Corollary 12. II is not valid in the case p = 1. To see it, consider the measure: 

j 

where Xj = Xj(t) G C l (I;£l) and ay G C]_(I,W) such that J2j a j( t ) = 1 Vi € i". We can 
approximate fi by a sequence of measures /i m G Ai(/io, fii) as follows: For each m G N consider 
the division i[ = /c/m, < k < m of /. Let be the optimal solution of Kantorovich 
problem due to W^^p^), and T® k := Id + (t - t< m) ) [T ro>fc - Id] /(t { ™\ - t™). 
Define fj, m as follows: 

Mm,(t fc ) - ^(tk) J Mm,(t) - 1 m ,fc,#/V(*k) > r fc — — fc+l ' 

Then, by Proposition I2.1| // m are bounded in Hi and — ► M- However, p, G" Hi unless 
are constants in t. To see it, note that the continuity equation takes the form 

= ^/ (uj(t)M%j(t),t)+v j {t)-V x (t>(x j (t),t)dt = ^ J -a j( f)(x j (t),t)^v j (t)-x j ]-V x (l>(x j (t) 
j 1 j 1 

where Vj(t) are the velocities attributed to Xj. It is evident that, unless ay = 0, for any 
possible choice of vj one can find (j) = <f)(x, t) for which the integral on the right does not 
vanish. 
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3 Main results 

Let the pressure P = P(x,t) € C 1 (f2/) and the associated action: 



== ^IHli+ / P^{dxdt) ; ^eH 2 . (3.1) 
Let us recall the definition of the action Jp: 

Jp(x,y,h,t 2 ) = inf 1^L + P(x(t),t) dt ; x : [ti, t 2 ] fi, s(ti) = x , x(t 2 ) = j/| 

(3-2) 

Remark: iVote i/iai Jp is not a function on = R n /Z n in eac/i o/ £/ie variables x, y, 
separately. However, for each q £ Z n and each x, y £ M. n , ti,t 2 S /, Jp(x + g, y + g, ii, i 2 ) = 
J P (x,y,t 1 ,t 2 ). 

Definition 3.1. (L) (the relaxed Lagrangian): 

C{fio,m) := inf Lp(fi) . 

Me A2 (MO, Ml) 

Definition 3.2. (M). (the Monge problem): 

M{no,m):= inf / J P (x, T(x), 0, T)n {dx) . 
T # M0=Mi Jn 

Definition 3.3. (K). (the Kantorovich problem): 

/C(// ,Mi) := hif / J P (x,y,0,T)X(dxdy) 
A J(] 

among all probability measures on Q x Q with the same f2 marginals fJ-o^i- 
We now introduce the Hamilton-Jacobi (HJ) equation 

f + i|V^ = P. (3.3) 
Let us denote the set of classical sub-solutions of the H-J equation as 

A P := j^eC 1 ^) ; & + ^|V x 0| 2 <PV(i,i) in fljj . 

For our purpose we need a generalization of the concept of a classical sub-solution. The 
concept of viscosity sub-solution (see, e.g. [E]) is too restrictive for us. So, we define a 
generalized sub-solution of the H-J equation as follows: 

The set of generalized sub solution of the H-J equation is given by 



T P := {<p £ L/P(J2 7 ) ; Vx G C\P, SI), 
d_ 
~dt 



<f>(x(t),t) <^|x"(t)| 2 + P(x(t),t) holds for Lebesgue a.e t € /J (3.4) 
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Remark (i): Note that (p(x(t),t) is a Lipschitz function on / if is Lipschitz and x G C 1 (I). 
Hence it is a.e. differentiable (as a function of t) on / by Rademacher's Theorem (see, e.g., 
[E])- 

Remark (ii): It is not difficult to see that any classical sub-solution is also generalized sub 
solution, so Ap C Ap. The concept of generalized sub-solution is more general than that of 
a viscosity sub-solution. The relation between generalized sub-solutions and viscosity (and 
anti- viscosity) sub-solutions is discussed in section 15.11 

Definition 3.4. (E): (The Euler formulation) : 

S(jjto, fix) := sup <^ / 4>(x,T)nx(dx) - / 4>[x,Q)fi Q (dx) 
4>eA*(P) 

We now state our main result: 
Main Theorem: 

Assume P G C^fij). For any /x , Ml G M.: 

£(Mo, Mi) = £(Mo, Mi) = £(mo, Mi) • (3-5) 

There exists minimizers \i G A2(/Z0)Mi) °/L (Definition Iff. i)) and a maximizer ifj G Ap o/E 
(Definition \3.4\) such that 

il>t + \\Vx4>\ 2 = P ;a.eonfi 7 . (3.6) 

Assume, in addition, there exists C(t) > on Iq so that P(x,t) — C(t)\x\ 2 is a concave 
function on W 1 for any t G Iq. Then, for maximizer ip of £(mo> Ml)? ^ ere exists a closed set 
K C 0/ suc/t i/iai 

i) TTie restriction of if) to Kq := KP\Vti is continuously differentiable, the equality (j3.6|) 

/joZds for any (x.t) G -Ko and V X V is Locally Lipschitz continuous on Kq. 

ii) Let v be a Lipschitz extension of \7 x ip to $7/ . Let T = be the flow generated by v. 

Then Kq is invariant under this flow. 

iii) A minimizer fi G A2(mo> Ml) °f L is n °t necessarily unique. However, any such minimizer 
is supported in K and the vectorfield v = V x ip is uniquely defined on the support of 
any such minimizer. 

iv) Any such minimizer is transported by the flow T, that is 

[ T t?] # M(ti) =M(t 2 ) 

holds for any tx,t 2 G (0,T). Moreover, if t x = (res. t 2 = T) then T£ := lim T ^ T* 
(res. Tf = lim r ^T'TJ~) are continuous maps transporting /jq to [iu\ (res. fim to fix). 

v) The map are optimal with respect to the cost function c(x,y) = Jp(x,y,tx,t2) and 

the measures M(*i)' M(t 2 )> where either tx G /, £2 G Jo or ti £ Iq, t 2 G I. 

vi) If P = then the optimal solution ^ of E (Definition 13.4(1 is in Cj^(Qj ). In particular, 
the flow T can be defined anywhere in terms of tp as 

T\\{x) = x + (t 2 -tx)V x 4>(x,tx) ,Vti < ta G I ,Vx G O . 
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4 Dual representation 

The key duality argument for minimizing convex functionals under affine constraints is sum- 
marized in the following proposition whose proof is given in the appendix: 

Proposition 4.1. Let C a real Banach space and C* the its dual. Denote the duality C = C* 
relation by < c*,c >€ R. Let Z a subspace of C and h G C*. Let Z* C C* given by the 
condition z* G Z* iff < z* — h, z >= for any z G Z. Let J 7 : C* -* R U {°°} a convex 
function and 

I : = inf .F(c*) . 
c *ez* 

Assume further that Aq := {c* G C* ; .F(c*) < /} zs compact (in the *— topology ofC*). 
Then 

sup inf [.F(c*)- < c*, 2 > + < h, z >} = I . 
zgZ c*ec* 

In particular, both sides equal oo if Z* = 0. 
4.1 Dual representation of H 2 

We shall apply Proposition 14. II were the space C is all the continuous functions q = q(x,t,v) 
on fli x R 71 subjected to: 

\\q\\ := sup <{ ' V / [ 1 2 } <oo . (4.1) 



(z,M)en 7 xiR n 



l + \v 



The dual space C* contains all finite Borel measures v on Qj x R n of finite second moments: 



|i/|(cfe<iicfo) < oo ; / |t>| |^|(dx(it(if) < oo . 
Define the subspaces Z, Zo of C as 

Z :={z = <j> t + v-V x <j> ; ^C 1 (O / )}cZ:={z = ^ + J )-V^ ; <f> G C 1 ^) n LJP(0/)} 
Given /io,Mi £ A4, define hn ui as a linear functional on Z as follows: 

<V>^i> z>:= / 4>{x,T)m(dx) - \ <j)(x,0)no(dx) for z£Z, (4.2) 

Jn 

(in particular, < fyu 0)W , 2 >= if z G Zo). 

Lemma 4.1. T/ie functional hn n i; so defined, is continuous (bounded) on C. 

Proof. Let A be a probability distribution on x so that 7rj? A = /xo, ?r£ A = /xj. Then 

/ <t){x,T)n x (dx) - / ^,0> (^)= / / [^(y,r)-^(x,0)]A(dxdj/) . (4.3) 
Jn J Jnxn 
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Now, for £(s) : = 



(T-s)x+sy 



we obtain 



0(y,T)-0(z,O) 



(j)(C(s),s)ds 



V — x „ 



C(*),- 



In particular, 



|0(y, T) — 0)| < max max \z(x,t,v)\ < \\z\ 

(x,t)eQi \v\<Diam(Cl/T) 



z[as),s,V-A)ds (4.4) 



1 + /Diam(J2)\ 2 



where we used the definition on the norm || • || on C given by (|4.1|) . The proof follows from 
(j4.2l4.3f) and since A is a probability distribution on x 0,. □ 



The corresponding dual spaces are given by 



Z* := \u £ C*; / z(x,t,v)v(dxdtdv) = ,Vz G Z i 



(4.5) 



3 Z^o^ C*; 



z(x,t,v)v(dxdtdv) =< /i MOiAtl ,z > ,Vz G Z > . 



For any € H2, a convex subset of C* is given by 



i/eC; / <f>(x,t)v(dx,dt,dv) = / <f>(x,t)fj,(dxdt) V G C(Q/) 

f2 7 xR n JO/ 



Finally, F„ : C* 



U {00} is defined by 

w = 



5 /n,xM- MM^tftdv) if ^ G C* 

if^c; t 



OC 



We obtain 

Lemma 4.2. T/ie function is convex on C* for any [i £ H2. In addition, if F^iy) < 00 
and v £ Zq then v is a lifting of [i. Similarly, if F^u) < 00 and v £ Z* „ x i/ien /i G 

A2(mo,A*i)- 

Proof. The proof of Lemma 14.21 is almost evident from the definitions. Let us just prove the 
last part. Since Z* _ C Zq it follows that ^ is a lifting of fi £ H2. We only have to show 
that n £ A 2 (/i ,/Ui).' Let (f> £ C^Q/J n LJP(Q/), r? = r/(t) G Cft(I) satisfies < 77 < 1 on I 
and, for some e > 0, r/(i) = 1 for e < t < T — e, and rjt > on [0, e], r/t < on [T — e, T]. Set 
0( £ ) = ^on!lx[e,T-£] and 0( £ ) (x, t) = <f>(x, e) on t £ [0, e] (res. (/)& (x, t) = c/)(x, T — e) on 
t G [T- e,T]). Then 77^ G Cg(nj), so 



(# (£) )t + r/v • Va.<^ (e ) v(dxdtdv) 



T-e 



[0f + ?j • V^] vu\{dxdv)dt 



14 



+ / / Vt4>{ x i £ )V(t)(dx)dt+ / / r] t (p(x,T - e)fin)(dx)dt 
Jo Jn JT-eJn 

rjv ■ Vx<P(x,e)uM(dxdv)dt + / / rjv ■ V x (j)(x,T — e)v^ > {dxdv)dt . (4.6) 
JUxR n Jt-eJuxr™ 

Since v is a lifting of some /x G H2 it follows that vt t \(dxdv) is a probability measure 
on 17 x M n . By the Cauchy-Schwartz inequality we estimate the last two integrals by 
2|| V a; ^>||oo\/E !/ (H 2 )£ 1 / 2 . By Lemma 12.21 pLt t \ is Holder continuous of exponent 1/2 in t, 
with respect to the W\ topology, so 

n7]t4>{x 1 e)n ( t ) (dx)dt = [ [ r} t 4){x 1 e)^ ( Q ) (dx)dt + 0(£ l/2 )\\V x (t)\\ O c f 1 7ft | eft 
a Jo Jn Jo 



= ( cP(x,e)n {0) (dx) + O(e 1/2 )||V ;E 0|| OO / m dt = [ <f>(x,e)n {0) (dx) + 0{s L ' 2 )\\V x ^\\ o , 
Jn Jo Jn 

(4-7) 

using rjt > on [0, e], hence J Q £ \rjt\ = f^Vt = 1- Similarly 

[ T [ (v4>)mt)(dx)dt = - [<Kx,T- s) m (dx) + 0(eV2) II V^Hoo. (4.8) 
Jt-e Jn Jn 

Letting e^Owe obtain from 14.714.5)) : 

[<pt + v ■ V x 4>] v{dxdtdv) — j [c6(x, T)ii( T }(dx) — (f)(x,0)fi^(dx)] = . 

The above is valid for any <fi G C (Oj ) n LIP(Clj). Since v G Z* ^ by assumption, it follows 
that /U(o) = Mo an d H(T) = /J-i, hence fj, G A2(y«o,/ii). □ 

Corollary 4.1. I//16 H2 i/ien 



-||m||2 = - inf / + |V s 0| 2 /2)//(<tedi) = sup V 7 ^ . (4.9) 

2 0eCg J 2 ^eC* Jn x I V^l^e&ccft) 

as we// as 



( f Qi (f>tn(dxdt) 



inf \ / (ct> t + \V x (l>\ 2 /2)fi{dxdt)+ / 0(x,O)/io(da;)- / e/>(x,7>i(efe) 



j 5IHI2 */ M G A 2 (/io,Mi) / 4 1Q x 

1 00 «'/ M £ A 2 (// ,Mi) 
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Proof. Certainly, satisfies all the conditions of Proposition 14.11 Using Lemma 14,21 and 
Proposition 14.11 in the definition of [|/x[|2 (Definition 12. II for p = 2) we obtain that 

dHll = i?f Ff,(u) = sup inf {F^{u)~ <u,z>) 



sup inf 



^H 2 - <j>t-v- V x (, 



sup 



sup inf 

^ + ^|V^| 2 



~\v- VA\ 2 -&-\\V x (t>\ 2 



u{dxdtdv) 

u{dxdtdv) 



[i(dxdt) H — inf / \v — V 
2 vec; Jn I xR n 



x cj)'\ i v{dxdtdv) 



So, we set v = fj,5 v ^\/ x( f > to annihilate the second integral and obtain the first equality in (|4,9j) . 
For the second equality in ()4.9|) we observe 

inf {/ (<p t + \V x <P\ 2 /2)f,(dxdt)\ = inf inf ( / + /3 2 |V^| 2 /2)/x(^)l 



inf 

0€Cg 



4>t^{dxdt) 



2 / nj |V^| 2 /x(^) 

Finally, we obtain (|4.10|) using the constraint Z* Q „ a for Z* in Proposition 14.11 



V 



□ 



Example: Let /i = ^f =1 Pk^(x-x k (t)) where x^ '■ I — > satisfies J T |xj| 2 dt := | 1 2 < 00 



and $ > 0, £\ /% = 1. Then 

/ 4> t fi(dxdt) =y2@j / 



0* 



(xj(t),t)dt 



and 

On the other hand 

fT 



[ \V x <j)\ 2 Li(dxdt) = V/? i [ {V^ixjit),^ 
Jn T ■ Jo 



t)dt 



— (xj(t),t) - ±j(t) ■ V x (j)(xj(t),t) 



dt 



Xj(t) ■ V x (f)(xj(t),t)dt 
so, by an application (twice) of the Cauchy-Schwartz inequality, 



(/n, <f>tK dxdt )) (Ej Pj I T ^x<t>{xj{t),t) ■ Xj dt 



fsij N x cf>\ 2 Kdxdt) ~ £ 3 Pi lo \V x <t>{x3{t),t)\Ht 



< 



\xj\ 2 dt 
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In fact, it can be shown that \\nW2 coincides with the above sum, and that there exists a 
maximizing sequence <j) n {x,t) such that V x 4> n (xj(t),t) — > Xj(t) for all j and a.e t G I (even if 
some of the orbits Xj intersect (!)-see [W]). 

4.2 Dual representation of the Lagrangian 

We shall now define a strong convex perturbation of the Lagrangian Lp fDefinition l3.1|) . Let 
also F : R -> R+ U {00} such that 

F(q) = 00 if q < ; F(0) =0 ; of < F(q) < Cq w if q > (4.11) 

where 1 < u < 1 + l/(n + 1) and c, C > 0. The functional if : C* — > R U {00} is defined by: 

%f{v):= j eF(f u )dxdtdv + - I \v\ 2 v{dxdtdv)+ I P(x,t)u(dxdtdv), (4.12) 

if 1/ = f u (x,t,v)dxdtdv is absolutely continuous with respect to Lebesgue measure and the 
density /„ satisfies F(/„) G L 1 (f2/ x R n ). Otherwise Tf(y) = 00. Note that, since F(g) = 00 
for q < 0, it follows that lf{y) = 00 if ^ G C* is not a non-negative measure. However, Zf 
can attain a finite value also for a measure 1/ which is not normalized (i.e not a probability 
measure on fi/ x R n ). 

Given yUOj/^l £ -M, define 

lf(/i ,Mi):= inf TfH. (4.13) 

" eZ M .Mi 

Next, we claim 

Lemma 4.3. For any e > 0, 

if (fJLo,Hl) > C(fJLo,(ll) 

where >C(/xo, as Definition ^. 1\ 

Proof. First, we can restrict ourselves to non-negative measures v G Z* „ x , since otherwise 
J F(f u ) = 00 by (|4.11|) . We only have to show that if 1^ > and f G Z* „ x then is a lifting 
of some /U G A2(//o>A i i)- 

Using Lemma l4.2l it is, therefore, enough to prove that vt{dxdv) is a probability measure 
on x R n for a.e (Borel) tel. Setting <j)(x,t) = rj(t) G Cq(I) we obtain from (|4~5|) that 

1(1 u t (dxdv)) %t = 
Ji \JnxR n J dt 

for any such 77. This implies that f^ xRn v t {dxdv) is constant for a.e. tel. Since v > it 
implies that z/< is a constant multiple of some probability measure on f2 x R n for a.e. tel. 
This constant equals one since the fi marginal of v% is C* continuous on I by Lemma 12.21 and 
is a probability measure at t = (/^o) and t = T (m). □ 
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We now proceed to a dual formulation of the constraint minimization of If. Certainly 
if satisfies the assumption on T introduced in Proposition 14.11 In fact, it follows that the 
set {v G C* ; lf{y) < C} is bounded (and hence *— compact) for any real C. Then 
Proposition 14. II and (|4,13|) yield 

if (Mo,Ml) = sup inf \lf(y)- < u, z > + < h^^z >] 



dxdtdv+ 



sup inf 



eF(f)-f[<p t + v-V x <p--\v\ I -P 



(f>(x,T)(j,i(dx) - 4>(x,0)ii (dx) , 
where infy stands for the infimum on all measurable functions on f2/ x M. n . Let 

1, 



H e {fA) ■-- 



dxdtdv 



eF{f) 



\v\~P f 



'fi/Xi?™ 

Let F* be the Legendre transform of F: 



eF(f)-[<p t + -\V x 



dxdtdv . 



P)f 



F*{\) = sup [sX - F(s)] . 



By our assumption we know that F* is also convex and non-negative on R. It satisfies 
F*(A) = for A < 0. Now, 



inftf £ (/, 



A+n/2 



!i,xl" 



h .. t *+w/*-p_& hlrillllr 



Let 



and 



9M := -e^ 2 f G (— 

Jfl! \ 



G(s) := / F*{s - \v\ 2 /2)dv 
+ \V x cj ) \ 2 /2-P 



(4.14) 



We have proved: 
Lemma 4.4. For e > and G -M, 

if {^0,^1 



dxdt + I 4>(x, T)ni{dx) — 4>(x, 0)^o(dx) 
n 

(4.15) 



sup tf e (0) . 

4>eC 1 (Q Io )nLlP(n I ) 
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We shall also need the following result, whose proof is direct and omitted: 

Lemma 4.5. If F satisfies \4-H\) then, for some constant c > 0, the function G defined in 
[fTfl ) satisfies cq"/^ 1 < G(q) < c~ 1 (f/^~ 1 . Thus, the first integral of UUfy is estimated 
by 



.l+n/2 



G 



b t + \V x <P\ 2 /2-P 



dxdt 



-0(eT 



+ 



|V^| 2 



p 



dxdt 



where a = 1/(uj - 1) - n/2 > and s = u/{u - 1) > 1 + n (c.f. 
We also need: 

Lemma 4.6. Let /io>Mi £ M.. Then there exists a connecting orbit p G A2(po,pi) of finite 
H2 norm and a lifting v such that both \x and v has densities in L P (J7/) (res. L P (J7/ x W 1 )), 
where 1 < p < 1 + 1/n. 

In particular, it follows that for such v as guaranteed in Lemma 14.61 each of the integrals 
in Q4.12JI is finite. Hence, there exists C > (independent of e) and v G Z* such that 
Tf{y) < C for any e > 0. In particular, lf{p,Q,p,{) < C for any such e by (|4.13j) . It follows 
from this, Lemma I4..SI and Lemma 14.41 that 

Corollary 4.2. For any po,Pi G M there exists C > independent of e where 

C > sup # E (» > C(jio,iii) . 

^ec 1 (n /o )ni7P(n / ) 

Lemma l4.Hl is a direct result from Lemma 14.71 below. For its presentation we define the 
space Hp([to,£i]) by restricting H p = H p (7) to orbits defined for a time interval [to, ii]. The 
norm of p G H p ([to,ii]) is denoted by HMlkjto.til- Lemma f4.7l is also used in the proof of 
Lemma 15.81 

Lemma 4.7. For any t% > t<i G any path x = x(t) : [to, ti] — > ft, and any a > there exists 
an orbit pmdt G H p ([to,£i]) with pr to \ = 5 X0 , pu-A = S X1 where x% = x(ti), i = 0, 1, such that 



Mll2.Ito.t1] ^ / M dt + C\ti-t \a 



and p^(dx) = p(x,t)dx where p G L p (f2 x [i ,ii]) for any p G [1, 1 + l/ n )- Moreover 



\P\ P < C( P ) 



(ti - to) 



a 



t\ — to 



n(p-l)' 



VP 



supp(p) C { (x,t) G n x [to,ti] ; - z| < C- — — V t G [t ,*i] 



and 



In particular, the choice x(t) = xq + (%i — Xq) yields 
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a 



Ml2,[*o,ti; 



< 



t\ — to 

The proof of Lemma 14.71 is given in the Appendix. 



+ C\t 1 -t \a 
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5 Proof of main results 



5.1 On the Hamilton- Jacobi Equation 

In this section we introduce some fundamental results for the H J equation 

4>t + \W x ^\ 2 =P (x,t)€fi/ (5.1) 

where P G C l {Qi). The book of L.Evans [E] contains a detailed exposition on the Hamilton- 
Jacobi equation. However, the discussion in [E] is restricted to generalized solutions of viscos- 
ity type and for time independent Hamiltonians, which excludes the application of backward 
solutions and time dependent pressure P = P(x,t) . The results in this section are all needed 
for the proof of the Main Theorem in section 

We list below some properties of the action Jp (|3.2j) : 

Lemma 5.1. For P G C l {yti), the action Jp is satisfies the following: 

(a) For n < T2 G [0,T] and xx,x 2 G M n , the value of the action Jp{x\, X2, t\, t%) is realized 

along a (possibly not unique) orbit x which satisfies the equation 

W = V x P(x(t),t) ; t G [n,T 2 ] . (5.2) 

(b) Assume further that there exists C(t) > so that P(x, t) — C(t)\x\ 2 is a concave function 

on 1" for any t G Iq. Let x be an optimizer orbit connecting xi,r\ to X2,T2. For any 
y GR" andte (0,T) 

Jp(x 1 ,X2,n,T 2 ) - J P (xi,y,Ti,t) > 

x(t 2 ) ■ (x 2 - y) + 



P(X2,T 2 ) - -|a;|(T2 



(t 2 -t)- 0(\x 2 - y\ 2 ) - o(t - T 2 ) , (5.3) 

Jp(y,X2,t,T 2 ) - Jp{xi,X2,Tl,T 2 ) < 



x(n) ■ (x 1 -y) + 



P( Xli r 1 )--\x\(r 1 ) 



(n-^ + Odxi-yl^+oit-n) . (5.4) 



(c) For any x\, y, X2 G R n , t\ < r < t 2 

J P (x±,y,ti,T) + Jp(y,x 2 ,r,t 2 ) > Jp{xi,X2,h,t 2 ) 

holds . 



(5.5) 



(d) For any pair x\,x 2 G W 1 and a triple t\ < r < t 2 there exists y* G M n (possibly 
non-unique) for which the equality holds in (|5.5|) : 



Jp(xi, y*, h, t) + Jp(y*,x 2 ,T, t 2 ) = Jp(xi,X2,h,t 2 ) 



(5.6) 



There exists a (possibly non-unique) optimal orbit x connecting (xi,ti) to (x 2 ,t2) such 
that x(t) = y. However, for any such optimal orbit, x(t) is determined uniquely. 
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(e) For any t > t\ x £E W 1 and a.e y £ 







1 



Jp(x,y,ti,t) + -\V y Jp(x,y,ti,t)\ = P(y,t) . 



(5.7) 



Definition 5.1. <f>(x,t) is a forward solution of (|5.1j) iff, for any i£(! and ti > to 6 I 

4>(x,t 1 )= inf [Jp(j/,ar,tb,ti) + 0(i/,to)] . (F) 

y€K n 

Likewise, 4> is a backward solution iff 

<p(x,to) = sup [-Jp(x,y,to,h) +<l>(y,ti)] . (B) 



Remark: It follows, by the remark proceeding (|3.2I) . that the right sides of (F) (res. (B)) 
defines a function which is Z n periodic on W 1 , namely defined on Q, if </>(•, to) (res. (/)(■, t\)) 
is a function on £1. 

For the special case of zero-pressure Hamilton-Jacobi equation, the action is reduced to 



Jo(xi,x 2 ,t 1 ,t 2 ) 



\X 2 - X\\ 



2(t 2 -*i) 

and definition 15.11 reduces to the (original) Hopf-Lax formula: 

Definition 5.2. A forward solution of the pressureless Hamilton-Jacobi equation 







satisfies, for any t\ > to and if!! 

Six, U) = inf 
while a backward solution satisfies 

4>(x,t ) = sup 



1 \x — yY 

2 t\ — to 



1 \x — y\ 



2 t± — to 



+ 4>(y,to) 



(Fo) 



(Bo) 



A forward (backward) solution can be constructed from an initial (end) data at t = 
(t = T) as follows. 

Lemma 5.2. For any continuous initial data 4>(,0) on f2 and P G LIP(Qi), 

<f>(x, t) = inf [J P (y, x, 0, t) + (f>(y, 0)] (5.8) 



(5.9) 



is a forward solution and satisfies ()5.1|) a.e. Moreover, <p 6 LIP(VL x (0,T]) and 

|^(x,t) < £ 

| a; — y\ ~ t 
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where C is a constant independent on (ft(,0). Likewise, for any continuous end data <ft(,T) 

4>{x,t) = sup [-J P {x,y,0,t) +4>(y,l)] 



is a backward solution and satisfies (|5,1|) a.e., (ft E LIP(Q x [0, T)) and 

\(ft(x,t)-(ft(y,t)\ < _C_ 

\x-y\ ~ T-t { ' ' 

If, in either cases, the end data (ft(,0) (res. (ft(,T)) is Lipschitz on £1, then the corresponding 
forward (backward) solution is in LIP{Qi). 

Next, we establish the connection between generalized sub-solutions, as defined in (|3.4|) . 
and forward/backward solutions: 

Lemma 5.3. Both forward and backward solutions are generalized sub-solutions in the sense 
of (|3.4[) . A forward (backward) solution is a maximal (minimal) generalized sub-solution 
in the following sense: If ip is a generalized sub-solution and (ft is a forward (backward) 
solution so that ip(x,r) = (ft(x,r) for all x E Q and some r E I, then (ft(x,t) > ip(x,t) 
(4>(x, t) < ift(x, t) ) for all x £ £1 and t > r (t < t) in I. 

An immediate corollary from Lemma 15.31 is: 

Corollary 5.1. Let (ft be a forward solution and ip is a backward solution on $1/. 

i) If (ft(x, T) = ifi(x, T) holds Vx G 0, then ift(x, t) < (ft(x, t) V(x, t) E Qi- 

ii) Similarly, if (ft(x,0) = ip(x,0) holds Vx E ft then ift{x,t) < (ft(x,t) V(x,t) E Jl/. 

Next, we wish to address the notion of a reversible solution: 

Definition 5.3. A reversible pair {(ft, (ft} where (ft ((ft) is a forward (backward) solution on 
such that (ft(x,0) = (ft(x,0) and (ft(x,T) = (ft(x,T) for any x E 0. By Corollaru \5.1\ (ft > (ft on 
0,1 . For any such reversible pair we denote the reversibility set of the pair as the relatively 
closed set Ko((ft,(ft) C Qi given by 

Ko((ft,(p) := {(x,t) E f2j ! Hx,t) = $(x,t)} C /o 

Likewise, 

K® (0, f) = K $, 0) n [n x {t}) for any t E (0, T) . 

If (ft = (ft then (ft := (ft = (ft is called a reversible solution. 

From Corollary 15.11 we obtain a way to create reversible pairs. It turns out that, in the 
case P = 0, this way yields reversible solutions: 

Lemma 5.4. Given (fto E LIP(Q), let (ft be the forward solution subjected to (ft(x,0) = (fto(x). 
Let ip be the backward solution subjected to ip(x,T) = (ft(x,T), and ift the forward solution 
subjected to 

Then {ip, ip} is a reversible pair. Moreover, if P = then ip = ip is a reversible solution. 
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The next Lemmas indicate that reversible pairs (in particular, reversible solutions) are 
closely related to classical solutions of the Hamilton-Jacobi equation. 

Lemma 5.5. If <fi £ C (Oj) is a classical solution of (|5.1[) then 4> is a reversible solution. 

Using Lemma 15. II we show that the converse of Lemma 15.51 also holds, in some sense: 

Lemma 5.6. If {4>, eft} is a reversible pair then both eft and <fi are differentiable on Kq := 
Ko((fr,4>) (cf, Dehnition \5. 6 J\) . Moreover, V0 := {V x ^, 4> t } = V<^> := {V x (/>, 4>} and the H.J 
equation is satisfied on this set. If, in addition, P satisfies the condition of Lemma \5.1Y (b) 
then V(p is locally Lipschitz continuous on Kq and <f> satisfies (|5.1|) pointwise on this set. 

Lemma 5.7. Assume P satisfies the condition of Lemma \5.1Y (b). Let v(x,t) be a Lipschitz 
extension ofV x (j) from Kq to fi/ . Then the set Kq is invariant with respect to the (unique) 
flow generated by the vectorfield v. 

Finally, we introduce the two following results, to be needed in Section T5.3I 

Lemma 5.8. If cj) G C 1 (r2 x [to^i]) satisfying 

<k + \\V x <i>\ 2 = P + Z ; (x,t)€fix [t ,ti] 

where P,£ € Lip(£l x [tQ, t\]), s > n + 1 and ||£|| s stands for the L s (f2 x [tQ, t\}) norm o/£ then, 
for any xq, x\ in Q, any t\ > tQ and any orbit x = x(t) : [to,ti] — * satisfying x{to) = xq, 
x(ti) = Xf. 

<t>[x x ,t x )-^xo,tQ) < \ r !^! 2 + f 1 ptm^dt+dMf (tx-to^+c^pwupih-tormf 



s > 

to J to 



where p - 2 p+n(p-i) > , P - s* .- —, A - 2p+np _ n , V - 2p +n(p-i) ■ 
From Lemma 15.81 and Definition 15. II we also obtain 

Corollary 5.2. Let <fi G be a solution and ip a forward solution of the respective 
equations on fij: 

<Pt + l/2\V x 4>\ 2 = P + i ■ ^ + 1/2|V^| 2 = P 
such that tfj(x,0) = (p(x,0) on 0. Then 

*P(x,T) > cf>(x,T) - [cM\\f + C 2 \\P\\ lip mf s 

where s, f3 as defined in Lemma f57%l 
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5.2 Proof of auxiliary results 

Proof, (of Lemma 15. 1(1 We shall only establish part (b), since all the rest follows by standard 
arguments. 

To establish this part, let us redefine the action Jp in terms of parameterized orbits 
x(s) — ► (y(s),r(s)), s £ [0, 1]. The optimal orbit is denoted x := {y, r}. The action takes the 
form 

1 



Jp(xi,x 2 ,t 1 ,t 2 ) = inf 

y< T Jo 



2r'(s) 



y(s) +P(y(s),T(s))r (s) 



ds 



(5.11) 



where .' stands for s derivative while {y, r} are the set of orbits satisfying y(0) = x\, y(l) = x 2 , 
r(0) = ti, t(1) = t 2 and r > on [0, 1] . An optimal orbit {y, r} satisfies the Euler-Lagrange 
equations: 



d_ (]_dy 

ds \t' ds 



d 


1 


y 


ds 


2 


i 

T 



P(y(«),r(s)) 



-P T (y,T)r 



(5.12) 



Now let us replace x 2 by y and t 2 by t. Deform the optimal orbit {y, r} into 

y(s) = y(s) + (y - x 2 )s ; f(s) =t(s) + (t - t 2 )s 

Now, {y,f} is an admissible orbit for the action with end-points (a?i,Ti) and (y, t). In 
addition, our assumption on P yields: 

P(jf(a), <?(*))?'(«) <P(y(s),T(«))^(*)+^'(«)^(y(«),r(s))(y-y) 

+ r'(s)P t (y(s), r(s))(f - r ) + (r - T>(y( g ), r(s)) + O 2 (|y - y|) + o(f - r) , (5.13) 
as well as 



\y 



I_'l2 _'/~' _'\ i_'|2/~' — '\ 

'y r , y (y - y ) |y I (t - r ) 



7- < + 



_ , +0^|y-x 2 |) + 0^-r 2 |), 

2 T 



2f' - 2r 

so, by substitution of {y,f}, integration by parts, (|5.11j) and (|5.12|) 

Jp{xi,y,n,t) < 





1 










Jo 


bf'(s) 





y(s) +P(y(s),f(s))f (a) 



Jp(xi, ac2,Ti ) T 2 )+ = 7 7 -T-y (l)-(y-x 2 )+ 

T 1) 



P(y(i),r(i)) 



2(r) 2 



(t-r 2 )+0 2 (|y-x 2 |)+ (|t-r 2 |) 



Using y(l) = x 2 , r(l) = r 2 and y /r = x we obtain (|5.3|) . The second inequality ()5.4|) is 
obtained similarly. 

□ 

Proof, (of lemma 15. 2(1 . 

Assuming given by 1(5.8(1 and ti > to ^ we have to prove (F) in Definition 15.11 namely 



<f>(y,h)= inf [Jp(w,y,t ,h) +(f>(w,t Q )] 



(5.14) 
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However, from (|5.8|) 

<t>{yM)= inf [Jp(z,y,0,ti) + <p(z,0)] and (j)(y,t ) = inf [Jp(z, y, 0, i ) + (f>(z, 0)] . 
from dHlI): 

Jp(z,y,0,ti) = inf {J P (z,-u;,0,to) + Jp(w,y,t ,ti)} 

u>GK n 

SO 

(j)(y,ti)= inf {Jp(z, w;,0, t ) + Jp{w,y,t ,ti) + 0(^,0)} 
= inf {0(w,i o ) + Jp(w,y,t ,t 1 )} 

and (|5.14j) follows. The proof for the backward equation is completely analogous. 

To prove the second part we proceed, as in the proof of Lemma 15. 11 to consider the orbit 

x(t) = x(t) + (y - x)r/t 

where x is an optimal orbit connecting (x,t) with (x(0),0), namely 



(x,t) 



\w\ /2 + P(x(r),r) dr + 0(5(0), 0) = J P (x(0), x,0,t) +<f>(x(0),0) . 



Since x(t) = y we obtain, using (|5,2j) and integration by parts 



/' 




/Q 





\x - y\ 
2t 



In particular we obtain that 



*fof) -*f*'*) < |^)| + + o(t 2 )) < £ 

\x-y\ ~ 1 Wl 2i v v - t 



{l+0(t 2 )) 



(5.15) 



where C = Tmax|x(i)| + Diam(Q) + 0(T 2 )]. Here max|x|(i) is the maximal possible 
value over all possible optimal orbits x. We now observe that there is a bound \x\ < D for 
any optimal orbit x where D depends only on the bounds of P and P x on Qj. This follows 
since there is a bound on the variation of x in terms of max(|P x |) via 1)5.2(1 . and a bound on 
the action Jp itself in terms of max(|P|) and max (|x|) by its definition. 

The reverse inequality of (|5.15|) and the result for the backward solutions follows analo- 
gously. □ 

Proof, (of Lemma 15 .HJI Let 4> be a forward solution. We know that is Lipschitz by 
Lemma 15.21 If (p is not a generalized sub-solution, then there exists e > and a C 1 or- 
bit x so that the set of points 

S t ; j t <j>(x{t),t) > M-+P(x(t),t)+s\ Cl 
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is of positive Lebesgue measure in I. Let tq be a density point of this set. Then, for 
t\ < tq < t% sufficiently close, and x\ = x{t\) ,x 2 = x(t 2 ), we obtain 



(x 2 ,t 2 ) - (f>(xi,ti) 



*2 ^ ['t2 

—4>(x(t),t)dt > 



\x(t)\ 2 



+ P(x(t),t) 



dt > Jp(xi,X2,h,t 2 ) 



which contradicts Definition 15.11 Hence, a forward solution is necessarily a generalized sub- 
solution. 

We now show that a forward solution <j) is a maximal generalized sub-solution. Let t/j be 
a generalized sub-solution, £ O/. Let r < t and x : [r, t] — > Q be a C 1 orbit satisfying 

x(i) = x. Assume (/> be a forward solution in 17 x [t, t] such that (f){x,r) = ip(x,r) on 17. 
Then, by definition 



r(.r./) :/-(.)'( r). -) + / -^-tp(x(s),s)ds < i/j(x(t),t) + 





1 1^ M 


r 









+ P(x(s),s)ds 



< 0(x(r), r) + J p (x(t),x, r, t) . (5.16) 

However, the same inequality survive if we take the infimum of the r.h.s of (|5.16|) over all 
such orbits x. By definition, this infimum is nothing but the value of the forward solution 
<f)(x,t). This implies the maximality of 4>. The statement for a backward solution is proved 
analogously. □ 

Proof, (of Lemma 15. 4 J) : 

Prom Corollary 15.11 and by construction it follows that tp < cj) and ip > i|) on SI;, while 
ip(x,0) = ip(x,0) for x S Q. On the other hand, since both ip and <p are forward solutions 
and 

^(x,0) =±(x,0) < <j) Q (x) 
on £1, then necessarily ip < <p on Qj. 

^(x,T) < 4>(x,T) = ip(x.T) 

and the equality ip(x,T) = ip(x,T) follows. 

Let us turn to the case P = 0. Set ipo = tp(-,0) (res. ipi = ip(-,T)). From (Fo), (Bo) in 
Definition O 



ipl{x) 
In particular 



inf 



\x - y\ 
2T 



tpo(x) = sup 



\x - y\ 
2T 



+ ipi(y) 



\x-y\ 
2T 



+ Mx)-Mv) > o 



(5.17) 



for any x, y £ R n . We may now lift the functions ipo,ipi from the torus O to periodic functions 
in M ra , and set 

2 2 

*o(x) = TMx) + y ; ^i(y) = y-T^i(y) , x,y£R n 
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as functions on M. n . Then (|5,17|) is equivalent to 

* (a?) + Vi(y) -x-y > ; Vx,i/Gl n . 



(5.18) 



The condition of reversible pair is manifested in being related by the Legendre trans- 

forms: 

V (x)=Bap[x-y-Vi(y)]:=Vt(x) ; ®i(x) = inf [x ■ y - tt (v)] := ■ 



In particular, both ^1 are convex functions on M. n . Recall that the super- gradient d x ^f of 
a function \& is given by 

y e d x ^ <=^ W(z) - *(ar) >y{z-x) ; Vz e M n 

and the equality in (|5,18|) holds iff y £ c^^o (equivalently, x E dy^i). 
Let now x,y,z £ M n , r S (0, T). The inequality 



I 12 I 1 9 i |2 

kc — w x — z\ \z — v\ 



2T ~ 2r 2(T-r) 
holds, and turns into an equality iff (T — t)x + ry = Tz. We obtain 



(5.19) 



tp(z,r) 



inf 



\z-y\ 

2t 



+ V'o(y) 



> t/>(z, t) = sup 



"2(T-r) 



+ V'i(y) 



^(z,t) - ip(z,r) 



inf 



I 1 9 I 1 2 

\z — x\ \z — y\ , , v , , s 



(5.20) 



2r ' 2(T-r) 

Let now l(z, t) := {x, y £ M. n ; (T — t)x + ry = Tz}. If we restrict the infimum in (|5.20|) to 
x,y £ l{z,r) then, by the equality in (|5.19|) . it is estimated by 



< inf 



\y - x\ 

2T 



+ i>o(x) - ij)x{y) 



T- 1 inf [^ {x) + ^ 1 {y)-x-y] 

x,y&(z,T) 



Now, the last term is zero iff there exists a pair (x, y) where both (T — r)x + ry = Tz and 
y £ d x ^o hold, namely: 



Tz — (T — t)x 



o 



-z := z T £ d x ^> 7 



T T 

where ^ T (x) := ^o( x ) + ^F" 3 - 2 - Finally, z T £ d x ^ T iff z T is in the domain of the Legendre 
transform Since is convex, it follows that \& T is strictly convex and it's Legendre 
transform is defined on W 1 . □ 

Proof, (of Lemma I5.5JI : 

Let x : [to, t±] — > R n be any orbit. Then 



/o 



Ixl 2 



+ P(x(t),t) 







dt= [ 




Jt 





+ ( t (x,t) + -|V^(x,t)| 
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'to 



-\T- V x (j)(x, t)\ 2 + (<f) t (x, t)+W- V x cp(x, t)) 



dt 



to 



±\ic-V x <p(x,t)\ 2 + ^ t <f)(x(t),t) 



1 f fl 

77 / \x - V x 4>{x,t)\ 2 dt + 4>(x 1 ,t l ) - <t>(x ,t ) 
1 J t 



where xq = x(to), x\ = x(t\). Hence 



to 



+ P(x(t),t) 



dt + (p(xo,t ) 



\x - V x (p{x,t)\ 2 dt 



to 



holds for any orbit x(t). In particular, 



0(aJl,*i) < 



to 



i^-+P(x(t),t) 



dt + 4>(x(t ),t Q ) . 



Moreover, if x is a solution of x = V ' x (j>(x, t), x(ti) = x±, to < t < ti, then the equality holds. 
Hence (j) is a forward solution. By the same way it follows that <j) is a backward solution, as 
well. □ 

Proof, (of Lemma 15 

Let r E (0, T) and Since is a forward solution and is a backward one, then 

(i)^{x,T)-0( yi ,O) < J P ( yi ,x,0,T) ; {ii) <j>(y 2 ,T)-4{x,T) < J P (x,y 2 ,r,T) V Vl ,y 2 eR n , 
Moreover, there exists y*,y 2 £ ^ n f° r which 

(i*)Hx,r) - ^,0) = J P (yl,x,0,T) ; (ii*) <j>{y* 2 , T) - 0(x, r) = J P (x, y 2 *, r, T) , 
while, for any x* G and i 6 (0, T) 

(i**) 0(x*,t)-^(i/J,O) < J P (yI,x*,0,t) ; (it**) <j>{y* 2 ,T) - $(x* ,t) < J P (x* ,y* 2 ,t,T) . 
From (**), (i**) and 

0(x,r) - 0(x*,t) > J P (y*,x,0,r) - J P (yJ, x*, 0, t) > 



xi(r) • (x - x*) + 



1 



P(x,r)--|x!(r)| 



{r _ t) _0 2 (\x-x*\)-o(\t-r\) (5.21) 



where xi : [0, r] — > R ra is an optimal orbit realizing the action Jp(y*, x, 0, r). Similarly, from 
(i**) («**) and JOJ 



(x, r) - 0(x*, t) < J P (x*,y 2 ,t, T) - J P (x, y 2 ,T, T) < 



x 2 (t) ■ (x - x*) + 



P(x,r)-l|x 2 (r)| 



(r-t) + 2 (\x-x*\) + o(\t-T\ 



(5.22) 



where x 2 : [t, T] — ► R n is an optimal orbit realizing the action J P (x,y 2 ,r,T). 
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Let now (x, r) G Kq 



Then Six, t) 



(x,t) = <P(x,t) and (l5~2T]) . (I5~22j) yield 



xi(t) • (x - X*) + 

x 2 (t) • (x - X*) + 
as well as 

xi(t) • (x - X*) + 

x 2 (t) • (x - X*) + 



P(x,t) - -|xi(r) 



{T _ t) _ \\x-x*\)-o{\T-t\) 



< 4>(x,t) - 4>(x*,t) < 



P(x,r)-i|x 2 (r)| 2 



P{ x ,r) - -^{r)] 2 



(T-t) + 2 (\x-x*\) + o(\T-t\) 



{T _ t) -0 2 (\x-X*\)-0(\t-T\) 



< 0(x,r) -0(x*,i) < 
1 



P(x,r)--|x 2 (r)| 2 



(r-t) + 0^(|x-x*|)+ (|t-r| 



Now, in order to show that both V x 0(x,r) and V x </>(x,t) exist for (x,r) G i£"o we only 
have to show that x"i(r) = x 2 (r). Adding (i) and (ii) we obtain: 

</>(y 2 ,T) - <t>( yi ,0) < Jp( yi ,x,0,T) + Jp(x,y 2 ,T,T) (5.23) 

for any yi,y 2 G W 1 . From (i*), (ii*) we obtain that equality holds in (|5.23f) for y*,y 2 : 

</>( y * 2 ,T)-<t>(y* 1 ,0) = J P (y* 1 ,x,0,T) + Jp(x,ylT,T) . (5.24) 

On the other hand, by pair-reversibility: 

<KvlT) - <f>(yt,0) < J P ( yi *,y 2 *,0,T) . 

Together with (f5~2l^) : 

Jp(y* 1 ,x,0,T) + J P (x,y* 2 ,T,T) < Jp(yl,y* 2 ,0,T) . (5.25) 

Comparing the above with (|5.5|) we obtain an equality in ()5.25|) , This implies that xi(r) = 
x 2 {t) by Lemma 15. ll -d. The existence of 



and of 



Vz0(x,t) = x(t) := xi(r) = x 2 (t) 



4t(x,T) =P(x,r)--|^| 2 (r) 



(5.26) 



follows from (|5.21|) . (|5.22|) . In particular the equality ()5.1|) holds for any (x,r) G Kq. 

We now prove the Lipschitz continuity of V x (/> on Kq. By Lemma 15.11 (b) we obtain 
that for any t G (0,T) there exists C = C(t) such that both Jp(y,x,0, t) — C(t)x 2 and 
Jp(x,y,t,T) — C(t)x 2 are concave functions of x for any y G M n . It follows by Definition l5.il 
(F) that q(x) := <ft(x, i)— C(i)x 2 , being the infimum of a set of concave functions Jp(y, x, 0, t)— 
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C(t)x 2 + cf>(y,0) ; y £ W 1 , is concave as well. Likewise, from (B) of the same definition 

we obtain that q(x) := 4>(x,t) + C(t)x 2 , being the supremum of a set of convex functions 
-J P (x, y, t, T) + C(t)x 2 + cp(x, T) ; y £ R n , is convex. 

Let x, x* £ Kq for some t £ Iq. Then both x, x* are differentiability points of the convex 
function q, hence 

(V x q-V x *q,x-x*)>0 . 
Using V<p = for x, x* we obtain 

{V x (p- V x *<f),x-x*) > -C\x-x*\ 2 . (5.27) 

The same argument applied to the concave function q yields 

(V x (j) - V x *<f>, x - x*) < C \x - x*\ 2 . (5.28) 

We may assume, by shifting and rotating the coordinate system, that x = (0, ... 0) and 
x* = (C, ... 0). Then (|5~2T1 KM yield 

|0,xi(0, ... 0, t) - ^(C, ... 0, *)| <C|C| • (5.29) 

It is left to be shown that there exists a constant C* such that 

|0 lX4 (O,O,...O,t)-^ a!i (C,O...O,t)| <C*|C| (5.30) 

holds for any i > 1. By subtracting an affine function from <p and 4> we may assume that (for 
the fixed value of t) 

^(O,t) = 0(O,t)=5(O) = g(O)=O (5.31) 

and 

V^(0, t) = V^(0, t) = V x q(0) = V x q{0) = , (5.32) 
hence ([5 .29(1 is reduced to 

|0 >a:i (C,O...O, t)| <C|C| , (5.33) 

and (|5,3U|) takes the form 

\MC,0,---0,t)\ <C*\C\ (5.34) 

for i > 1. 

Since q (q) is concave (convex), then 1)5.311 15~32|) implies that 

q < < q (5.35) 

everywhere. 

To proceed, we use the concavity of q to obtain 

q{(,y...0)-q(C,0...0)<q }X2 (C,0...0)y = (l )>X2 (C,0...0,t)y . (5.36) 
From ()5.36)) . ()5.35)) . the inequality <p > <j) and the relation between q and (f> we obtain 
<t>^ (C, ... 0, t)y > q(C, y, . . . 0) = 0(£ y, ... 0) - C(C 2 + y 2 ) 
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> <MC, v, • • • 0) - C(C 2 + y 2 ) = q(C, y, . . . 0) - 2C(C 2 + y 2 ) > -2C(( 2 + y 2 ) , 



so 



^ X2 ((,0...0,t)y + 2C\y\ 2 > -2CQ 2 
holds for any i/6R. Take now the minimum of the lhs above with respect to y to obtain 

X2 



^ > -2CC 2 



. jiE2 (C, 0, . . . 0, t)| <4C|C| 

and (|5.3(Jf) holds with C* = 4C where i = 2. The general case of i > 1 is evident. □ 
Proof, (of Lemma l5.7j) : Let (</>,</>) a reversible pair. We first note that, for any (x,t) G f2/ , 

< t) - #c, t) < 0) - T) + J P {y, x, 0, t) + J P {x, z, t, T) . 
However, from definition of Jp 

J P (y,x,0,t) + J P (x,z,t,T) < / - |5;(s)| 2 + P(x(s),s) ds , 

Jo l z 

holds for any C 1 orbit x which satisfies x(0) = y, x(t) = x and x(T) = z. It follows that 

1 2 „/_,., 



o < 4>(x,t) - 4>(x,t) < 





1 \H \\ 


f 


— \x(s)\ 


10 


2 1 ^1 



ds + (j)(x{0),0) - (j)(x{T),T) 



holds for any (x,t) G fi/ . 

As a result we conclude that, if there exists an orbit x crossing the point (x, t) G Qi such 
that 



^m S )\ 2 +p(x(s), S ) 



ds + 0(x(O), 0) - 4>(x(T),T) = , 



(5.37) 



then (j)(x,t) = (j)(x,t), i.e. (x,i) G i^o- Moreover, in that case, the entire orbit is contained in 
-Kcb namely U re / (x(r), r) C -Ko- I n addition, it must be a minimal orbit for the action and, 
by Lemma l5~Tl and Lemma EH V x 0(x(i),t) = V x (j)(x(t),t) = x(t). 

Next, let (x,t) G Kq and let x\ : [0,i] — > f2 (x~2 ■ [t,T] — > $7) the minimal orbits for 

which 4>{x,t) = </>(^i(0),0) + J * \ |xi(s)| 2 + P(xi(s), s) ds (res. 0(x,i) = (p(x2(T),T) — 

T ' 1 I * 1 2 ^^^^^ 

.ft 2 F 2 ( s )| + -P(^2(s), s) ds) and xi(t) = X2(t) = x. Then x = x\ U X2 satisfies (|5.37|) . In 
particular, it is contained in Kq and is an integral curve of the vectorfield V^. This implies 
the invariance of Kq under the flow of V x (j). □ 

5.3 Proof of the main Theorem 

First, the existence of a minimizer for £(/xo,/ii) in A2(/Uo,/xi) follows immediately by the 
lower-semi-continuity of ||^||2 and the local compactness of H2. Next, we shall prove the 
chain of inequalities: 

S(hq,h{) > £(/x ,/ii) > K(^Q,m) > £(/i ,/ii) 
from left to right, together with the existence of a maximizer for £(/zo,/ii) in A* P . 
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£{vo,Vi) > -C(/Uo,Mi) 

From Lemma 14.51 and Corollary 14.21 there exists a sequence e k — > and G 
such that 

£<* £fc (0 fc ) <-0(^ a ) H6II2+ / </>k{x,T)in(dx) - [ <f> k (x,0) fM) (dx)<C (5.38) 

Jn Jn 

where a>0, s>n + l and 

Zk = <i>k,t + \v x <l>k\ 2 /2-P ■ 

Let now S : O; ^ be a flow such that 



rT 

sw Pxen Jo 



dZ(x,t) 
dt 



dt := E < oo, 



ii) E(x, 0) = x, S # (-, T)hq = m 
By Lemma 15.81 and (i) we obtain 

fa(E(x, T),T) - tf» k (x, 0) < |e + |P|oo + Ci||e*||f + C 2 ||P|| % ||e 
Integrate the above against /j,q on 17 and use (ii) to obtain 

( (f>k(x,T)fj,i(dx) — ( <t> k (x,0)no(dx) < ^+|P| 00 + Ci||e fc ||f + C 2 ||P||i ip ||^||f 
isi Jn 1 

(5.39) 

Using s > n + 1 and /3 < 1/2 (c.f. Lemma l5.8|) we obtain from Q5.38JI and ()5.39|) that 
||6c||s — ► as e k — > 0. In addition 



lim inf 

A;— >oo 



4> k (x,T)/j, 1 (dx) - / fc (:r,O)^o(dx) 



> £ 



(5.40) 



Now, we may replace the sequence by a sequence of forward solutions ip k of the 
equation 

i>U + o l V ^fc| 2 = P ! V^O^O) = (j) k {x,0) . 



This is also a maximizing sequence which, by Corollary 15.21 together with ||^| 
yields 



lim inf 

k— too 



ip k (x 



,T)[ii(dx) - 

Jn 



i/j k (x,0)na(dx) 



> C . 



0. 



(5.41) 



From Lemma 15.21 we also obtain a uniform estimate on ifj k in LIP(Q x (0,T]). In 
particular, the sequence ip k (]T) is uniformly Lipschitz on f2. 

Now, define ip to be the backward solutions of ()3.6|) subjected to ip k (x,T) = ip k (x,T). 
From the first part of Lemma I5.21 ip (x, 0) < ip k (x,0) on $7 so (j5.41|) is satisfied for 
ip k as well. Moreover, by the last part of Lemma 15.21 tp are uniformly bounded in 
the Lipschitz norm on f2/. So, we can extract a subsequence of ip, which converges 
uniformly on LIP(£l x [0, T]) to a backward solution ^.In particular, both ip(,T) and 
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ip(,0) are Lipschitz. Let ip be the forward solution satisfying ^(,0) = ip(,0). By 
Corollary 15.11 and definition 15.31 the pair (tp,ip) is a reversible pair and both functions 
are in A P (see the first part of Lemma l5.3|) . Moreover, the inequality Q5.41JI is preserved 
in the limit process, so 



/ ip(x, T)p,\{dx) — / ip(x, 0)po(dx) > C 
Jn Jn 



holds for both ip = ip and ip = ip (recall ip = ip on Q x {0} and £1 x {T}). 
C(fio,(ii) > K{hq,ii\) 

Recall that there exists a minimizer of C(/iQ,pi) by the first part of the Theorem. Let 
p be such a minimizer. We now use the regularization Lemma 12.31 to approximate p by 
smooth densities p n = p n (x,t)dxdt. Let v n be the regularized velocity field. Then 



lim / \v n \ 2 (x,t)p n {x,t)dxdt = \\p\\\ 



(5.42) 



as well as 



Define 



lim / p n (x,t)P(x,t)dxdt = / Pp(dxdt) . 

JCt! J fit 



m n (x,t) = p n (x,t)v n (x,t) . 
Then m n £ C l {£li). Define now 

Pn{x,t) +e 

By assumption, v e n is Lipschitz on $7/, tel. Define pn\x,t) as the solution of 
dp^ 



(5.43) 



dt 



v nPn 



; p^(x,0)=p n (x,0) . 



(5.44) 



Since v e n is Lipschitz, we may define the flow associated with it as Tt* : 0, —* Q for 



t e I, namely T^{x) = y( x )(t) where y {x) = v e n (y^ x) (t),t) and y (x) (0) = x. It follows 
that r* e -j ^p n (-,0)dx = pn\-,t)dx for all tel. In particular: 

K (pn{x,0)dx,p&(x,T)dx} < J^ Pn (x,0)Jp (x,rj } (x)0,r) dx 
p n (x,Q)^-Jp (x,T*(x),0,t) dtdx 



p n {x,0) —Jp{x,T t {£) (x),t)+V y=n(x) Jp(x,y,0,t)-v £ n (r t £ (x),t) 



dxdt (5.45) 
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From If577|) 



d t J P {x, y, 0, t)+V y J P (x, y, 0, t)-v e n (y, t) = P{y, t)+- \v £ n \ 2 (y, t)-~ \V y J P (x, y, 0, t) - v £ n (y, t) 
Substitute the above in (|5.45l) at y = T*(x) to obtain 

tc(p n (x,0)dx,p^(x,T)dx) < p n (x,0) P{Tt(x),t) + i|<| 2 (r*(x),t) 



dxdt 



-\v^\ 2 (x,t)+P(x,t) 



dxdt < I p^Hx,t) 



~K| 2 (x,t) + P(x,t) 



where the last inequality follows from ()5.43j) . We next show that 



dxdt 
(5.46) 



lim 

£^0 



Pn( x >t) ~ Pn(x,t) 



dx = 



(5.47) 



for any t G I. In fact, we note that p n +e solves equation (|5.44|) . hence Wn^ '■= p n —Pn^+£ 
solves this equation as well. Since Wn\x,0) = e > we obtain that wffl > over f2/ 
and, moreover, 



Pn{x,t) - p%>(x,t) 



dx — \Q\e < 



w£\x,t) = [ ^i £) (x,0) 
Jn 



Me 



for all t £ I. Now we take first the limit e — > then the limit n — > oo in (|5,46|) . The r.h.s 
of (|5.46|) converges to C(po,pi). Now, p n (x,0)dx and p e n {x,T}dx converges, as n — > oo 
and e — ► 0, weak—* to po and pi, respectively. Since K, is lower-semi-continuous in 
both po and p\ we obtain the desired result from (|5.46f) . 



£{po,pi) < IC(po,pi). 



Let A E C*(Q x ri) be an optimizer of /C. Since 7r^ A = then 

(pi(x)pi(dx) = <f>i(y)\(dxdy) and / (f> (x) p (dx) = (f> (x)\(dxdy) 

n Jn Jn Jn Jn Jn 

for any continuous 0i, 02- Set 0i( 2; ) = ^(k^T) and (fto( x ) = VK^O) with ^ an optimal 
backward solution of problem £. Then 



1>{x,T)iJL X (dx) - [ 4>(x,0)p (dx) = [ [ [ij(y,T)-i>(x,0)]\(dxdy) . 
Jn Jn Jn 



Since ip is a backward solution then 



h%, r) - V(x, 0)] \(dxdy) < / / J P (x, y, 0, T)X{dxdy) = K . 
in Jn Jn Jn 
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We have proved 



i>(x,T)/ii(dx) - / tp(x,0)^ (dx) = / ■ip(x,T)fi 1 (dx) - / f_(x, 0)no(dx) = C(fi , Mi) • 
in Jn Jn Jn 

(5.48) 

We now turn to the proof of parts (i)-(vi) of the Theorem. 

i) Let fi^ be a minimizer of C. Given t G Iq, let //i/ 2 := M(i) ^ Let us consider 

to be the restriction of to x [0, t] and // 2 ) the restriction of /i^ - 1 to x [£, T]. 
Evidently, /Z 1 ) is a minimizer of Lp on the set of orbits A2(mq> M1/2) confined to [0,t] 
while /i*- 2 -* is a minimizer on A2(/ii/2, Mi) with respect to the same set, confined to [t, T]. 
In particular, 

Lp(^) + L P {^) = L P {^) = C{^^ 1 ) . (5.49) 
By what we know so far, 

^(x,t)fi 1/2 (dx) - [ ^(x,0)v o (dx) < Lp(^) (5.50) 
n Jn 

^(x,T)in(dx) - [ Jj(x,t)fi 1/2 (dx) < Lp(^) . (5.51) 
>n Jn 

However, if we sum (j5.50|) and (|5.51|) and use (|5.49f) and (|5.48|) . we conclude that there 
is, in fact, an equality in both Q5.5UJI and (|5.51|) . Same argument holds for ip as well. 
Thus ~~ 

i(}(x,t)nx/2(dx)- ip{x,0)no(dx) = Lp(/i (1) ) = / ijj_(x,t)fi 1 /2{dx)- ■tp(x,0)fj,o(dx) 
Jn Jn Jn 

Since ip(x,0) = ip(x,0), 

/ \^(x,t)-ip(x,t)]fj,u 2 (dx) = 0. 
Jn 

But, ip > tp by Lemma [5.31 Hence ip(x,t) = tp(x,t) on supp (M(t)J = supp ^1/2)- 
This, together with Lemma [5.61 proves that Supp (m^) H fi/ C Kq = {(x,t) ; t G 
Jo ,ifi(x,t) = ip(x,t)} and, in particular, that is differentiable at any point on the 
support of fi^ in SI/. 

ii) This part follows from Lemma 15.71 In addition, the limits lim T ^*r TJ and lim T ^o T^. 

exists since V x ip is uniformly bounded on Kq. In particular, the Lipschitz extension v 
can be chosen to be a uniformly bounded function on Qj as well. 

hi) Suppose there are two optimal solutions ipi,ip2 of £(mo ; Mi)- To prove the uniqueness for 
the vector field v we claim that 



/ IV^Vi-V^I 2 v(dxdt) = Q 
Jn! 
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for any minimizer fi £ A2(/xq, Mi) of C(/jlo, //i). Let ip = cr0i + (l — a)V>2 where a € (0, 1). 
Then 

IV^I 2 = a|V x Vi| 2 + (1 - a)|V^ 2 | 2 - a(l - a)|V^i - V^l 2 , 
so ^ + |V :c V| 2 /2 < P and 

/ + \V x ^\ 2 /2] fi{dxdt) < / P^(dxdt) (5.52) 
./nj " inj 

if X7 x ipi V x ip2 at some point in the support of a minimizer /j, (recall that both Vf/>j, 
i = 1,2 are continuous on the support of V> by Lemma l5~^)l . 

On the other hand, 

£(im,H\)= \ [ip(x,T)m(dx)-i/)(x,0)iJ,o(dx)]=£(no,fJ'i) (5-53) 

in 

follows from the assumptions that both ipi,ip2 are maximizers of £. Prom (J4,1U|) . (|5.52() 
and (|5.53|) it follows that 

M/-0 : = I P^(dxdt) > / ^ t +\^x^\ 2 /2)fi(dxdt) > £(no,Hi) = £{(*o 

* in/ 1 in 7 

in contradiction to the assumption that \i is a minimizer of -Lp. 

Let, again, fi S A2(/Uo,/Ui) a minimizer of C and 7/) a maximizer of £. Since ip satisfies 
the HJ equation on a closed set Kq containing the support of /u in Qj Q and is a C 1 
function there, we can extend it as a C 1 function on Qj so ip £ C 1 (J7/ ) Pi LIP(Qj) 
and, by (|535Jl . 



fi(dxdt)+ / [V>(a;,r)/ii((ix) - ^(x, 0)/j,o(dx)\ =£(^o,/^i) = -Lp 
in 



(5.54) 

We now use Corollary 14.11 Q4.1UJI to observe that ip is a maximizer of the left of Q5.54JI , 
so by taking the variation cj) = ip + er\ with r\ £ C (Q/) we obtain 

(jit + ^xtp ■ V x r))/j,(dxdt) + / ri(x,0)no(dx) - / i](x,T)ni(dx) > 

in in 

for any such n. Replacing r\ by —77 we obtain the equality above. Moreover, by the 
same argument following (|5.49l) to (|5.51|) we also obtain that 

/ (m + V x ip ■ V x rj)nt{dx) + / T](x,t )fj,, to ^(dx) - / rj(x, t)fj,^ (dx) = , (5.55) 
/t in in in 

hold for any < to < t < T. In particular, \i solves the weak form of the continuity 
equation with v = V x ^. 

Now, we know that, by the additional assumption on P, that Kq is invariant with 
respect to the flow Ti induced by the Lipschitz vectorfield v extending V^. We shall 
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now prove that fi is transported by this flow. That is, for any choice of to,t G (0,T), 
we need to show that fj,M = jt where 

It := [Ty #j Uto • 

Since T is the flow generated by v and Kq is invariant with respect to v it follows that 
7 = ~ftdt is supported on K and solves the weak form of the continuity equation as well. 
Setting Q = fi^ — j t , ( := Qdt we obtain from (|5.55j) 

/ / (rjr + v ■ V x ri)Cr{dx)dr= [ r](x,t)( t (dx) (5.56) 
Jt Jo, Jn 

for any r] E C 1 ([to,i];K) where we used Ct = 0. 

Let now h = h(x) G C 1 (f2). Let t] = 7](x, t) be a solution of 

r/ T + v ■ V x r] = ; r](x, t) = h(x) , to < t < t . (5.57) 

Since, by Lemma 15.61 and Lemma 15.71 the vector field V x ip = v is locally Lipschitz 
continuous on Kq which is invariant with respect to the induced flow, we can find a 
solution of 1(5371) onKn(!lx [to,t]) via 

r,(x,T) = h(T t T (x)) . (5.58) 

The function rj so defined can be extended into a C 1 function on $7 x [£o>£]- It satisfies 
(|5.57|) on Kq, so, recalling that ( is supported on Kq, we substitute now (|5.57|) in (|5.56|) 
to obtain Q = and the proof of part (iv) . 

v) The optimality of T is evident from the proof of (iii) and (iv). 

vi) From the last part of Lemma l5.4l it follows that tp is a reversible solution so Lemma 15.61 
implies that ipt -\-\S7 x il)\ 2 /2 = is satisfied everywhere on Qj . The flow induced by such 
a solution is given by T^(x) = x + (t — T)V x tp(x, r) and, by (iv) and (v), it transports 
^( T ) to optimally. 

6 Appendix 

Proof, of Proposition 14.11 
Define 

$(c*, z) := F{c*)- <c*,z> + <h,z> . 

First, note that 

/ = inf sup <3?(c*, z) . 

c*GC* zeZ 

Indeed, if c* Z* then sup 2eZ $(c*, z) = oo while, if c* G Z* then (f>(c*,z) = T{c*) by 
definition. We have, therefore, to show 

inf sup ^(c*, z) = sup inf Q(c*,z) . 
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It is trivial that 



inf sup 3>(c*, z) > sup inf <3?(c*, z) := I , 



inf sup $(c*,2:) < I . (6.1) 

c*GC* zg z 



so we only have to show that 
Define, for any z £ Z 

A 2 = {c* G C* ; *(c*,z) < /} . 
Note that (|6.1|) follows provided 

n a * + • ( 6 - 2 ) 

zGZ 

The next step is to show that, for any finite set z\, . . . z n £ Z, the set f] z . A Zi ^ 0. The proof 
of this part can be taken from the proof of Theorem 2.8.1 in [Ba]. 

Finally, note that Aq C Aq as defined in the Proposition, since I < I. It follows that Aq 
is compact, and that the non-empty intersection of finite sets implies (|6.2j) , □ 



Proof, of Lemma 14.71 

Let pi(r) be a smooth, positive function with compact support such that 



jn— 1 1 



r n l pi(r)dr = 1 ; 



r k pi{r)dr := M k ; 



-n-l 



p\dr := L(p) . 



Set also 
Define 



p a (r) = a n pi(ar) . 



v(x,t) 



?L^+W(t) ifto<t<(t + ti)/2 
^0.+W(t) if (to+ti)/2<t<ti 



p(x,t) 



(t-ta) 
1 



n Pa 



\x—x(t)\ 
t-t 



if to < t < (t + *i)/2 
if (to + h)/2 < t < h 



A direct calculation shows that p satisfies the weak form of the continuity equation: 

Pt + V x • (vp) = . 

Let us now consider the interval [to, (to + ti)/2]. The second interval [(to + ti)/2,ti] can be 
treated analogously. Define the lifting of p as 



/(x,t,t;)= ( T-^- 2 /" 



cxp 



(7 



p(x,t) . 



It follows immediately that 

/ v 2 f(x,t,v)dv = ^—p(x,t) + v 2 p(x,t) ; / |_f| p <fo=jf 



n/2^0.-p)/n fT n{X-p) 



a 



VOM) 
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Moreover: 

p{x,t)dx = 1 ; / pP(x,t)dx = (t-t ) n(1 " p) |S n_1 | / r n_1 pP(r)dr = a n ( p - 1 )(t-t ) n(1_p) |S n - 1 |L(p) 



|v(x,t)| 2 p(x,i)cte = \w(t)\ 2 + IS"" 1 ] / r" + Va(r)rfr = \w(t)\ 2 + IS^V^n+i 
In particular: 



"(t +*i)/2 /"(*o+tl)/2 



and 

r-(t +ti)/2 



v\ 2 f= / |^| dt + Ci|ti -t |a" 2 + 0(<t) 

/ = c 3 \ti - t \ n ^ +1 a n ^ 



in J t 

□ 



Proof, of Lemma 

We use Corollary 14, II with \x supported on $7 x [ii,io] and pt = 5 XQ , IHi = feu to obtain 



0(xi,ti) - 4>(x ,t ) < 



n J t 



(<t>t + \V x 4>\ 2 /2)pt {t {dx)dt 



By Lemma 14.71 we can find such a p for which 



< 



f l-|2 

/ a? | +Ci\ti — t \a 
J to 



and, for the density p = p^: 



P P <C 3 \t 1 -t \ 



n(l-p)+l a n(p-l) 



(6.3) 



(6.4) 



(6.5) 



n J to 



where p < 1 + 1/n and a any positive constant. Since p is supported, for any t, in a domain 
of diameter (ti — io)« _1 it follows 



n Jt 



+ |V^| 2 /2)/xfa;cft = O 



L(ti - t ) 5 



+ 



P(x(t),t)dt + 



to 



£pdxdt , 



n J t 



where L is the Lipschitz norm of P. By (|6.5j) we obtain 



^pdxdt 

in J t Q 

Collecting JOJ) to dnUl) 



< 



P liei| s < cl /p \h - t |Wi-P)+i]/p a -( P -D/P|| e || t 



(6.6) 



<Kx 1} tx)-<l>(x ,to) < I ! 1 l^l^t+C^IUti-tol^-^/V^-^+Cilti-tola-'+O (M^J^ 



The choice a = 1 1£| | s - 1 ) 7 where 7 - 2p+n(p _ 1) 
and yields the desired result. 



(w+lXg-l) 



and /3 



2p+n(p-l) 



is the optimal choice 
□ 
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